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i Abstract 

The Hubbard model in the limit of infinite U is investigated within a projected slave fermion representation and following a previous 
work of the author and collaborators [1]. The stability of the Nagaoka's phase with respect to a non vanishing concentration 
of holes (Sh) is analyzed by envisaging the existence of a spin effective action for itinerant magnetism of the Hubbard model. 
It is considered that, as the hole doping increases away from the half filled insulating limit, free holes are expected to be more 
delocalized. Depending on treatment for the hopping: a ferromagnetic or anti-ferromagnetic ordering might arise and the Nagaoka 
phase might have some stability with respect to 5h 7^ 0. 
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1. Introduction 

One of the few exact results for the Hubbard model (HM) 
is the ferromagnetic Nagaoka limit for U = oo |2l3j . This 
phase appears when one hole hops in the half filled band. 
Although this limit is not found in any material, it can 
be reached in optical traps [4j. The investigation of the 
Nagaoka mechanism provides relevant information about 
the phase diagram of the model and one starting point 
for understanding further the role of strong electronic cor- 
relations [1151617] and of realistic mechanism for itinerant 
(ferro)magnetism 5 3. In the limit of infinite Coulomb re- 
oulsion the HM is written as: 



the appearance and for the stability of the Nagaoka's phase 
following the long-wavelength approach with slave fermion 
representation worked out in Ref. pQ. By envisaging the 
derivation of a spin-effective action (S e ft) for the HM with 
very large U, a previous analysis was performed in Ref. pQ . 
On the other hand, in the present work, the role of delocal- 
ization of free holes close to half filling is investigated. It is 
considered that the increase of the number or concentra- 
tion of holes should increase the mobility of holes depart- 
ing from the half filled limit. An itinerant magnetic phase 
(ferromagnetic or anti-ferromagnetic) emerges depending 
on the structure and treatment of the hopping of spinless 
holes. 
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where Uj is a symmetric matrix with elements representing 
the hopping amplitude t only non-zero between the nearest- 
neighbor sites; Ci a is the projected electronic operator [5]. 
In this expression the chemical potential (i is to control 
the number of vacancies (away from half filling), and the 
projected electronic operator carries the effect of the strong 
correlations, i.e. it excludes the doubly occupied states. 

It is difficult to handle the strong electronic correlations 
and thus to provide exact results, in particular concerning 
the stability of the Nagaoka phase. However, a common 
trend is that this FM phase is unstable with respect to a 
finite concentration of holes in particular in the thermo- 
dynamic limit [8I9I10I11I12I13I14] . In the present work we 
investigate the role of the delocalization of (free) holes for 



2. Slave fermion for the U = oo Hubbard model 

To account for the strong correlations that forbid doubly 
occupied states, consider the slave-fermion decomposition 
for projected electronic operators [16] given by: cj CT = b\ a fi, 
where two operators have been used: b{ a stands for a spinon 
boson and f\ creates a charged spinless fermion. The func- 
tional generator for the U — oo Hubbard model is given by: 

Z v =oo =JV[b, b f ; /, / f ] exp (S^ [h,bl; f u //]) . The ac- 
tion can be written, with the time-dependent phase, as pQ: 
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The local non doubly occupancy (NDO) constraint is im- 
posed by a (local) Lagrange multiplier, Ai, by adding the 

term: H constr = \i(f}fi + J2a=t,l b L b ^ ~ 1)- 

The slave-fermion decomposition is equivalent to a par- 
ticular (lowest weight) representation of the su(2|l) super- 
symmetric projected electronic operators |6I17) : spinless 
holes are super-partners of spinons. A mapping for the vari- 
ables, incorporating implicitly the non-doubly-occupancy 
(NDO) constraint, is given by pQ: 
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and the corresponding variables for b\ and /• . With these 
new variables (zi,£i,fa respectively for bosonic spinons, 
spinless fermions and a local phase) the Lagrange multi- 
plier Aj is eliminated naturally, being the NDO constraint 
incorporated. With the spinon variables Zi, the images of 
the spin su{2) algebra - S - can be rewritten [1117] . for ex- 
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ample: Sf = g 1+ Lp ■ There is a local gauge invariance as 
consequence of the redundancy in parameterization of the 
electron operator in terms of the auxiliary boson/fermion 
fields. After some manipulation, decomposing the measure 
of the path integral into the new variables T>[b, tf; /, ft] — > 
Duspin (z,z)x D/j, fermion (£, 5 with the corresponding Ja- 
cobian [17] , the action reads [17|lj : 











S = y^ / iai(r)dT - ^ / & (<9 T + /x + Mi) ^dr 



' E' :i .' i "' < Zl ' z J > + hc ) dT 



(3) 



The first term of this action is a kinematical term and the 
second is the classical image of the Hamiltonian. The spin 
"kinetic" term (Berry phase) : ia = — < z\d t \z >= | fqprna i 
with |z) being the su(2) coherent state [1711] . The inner 
product of the su(2) coherent states is written 
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3. Factorization of the hopping 

By introducing more holes in the half filled HM, it can be 
expected they become progressively more delocalized. Con- 
sider that the band structure is such that the hopping term 
can be decomposed into two parts. One of them endows 
the holes with a dispersion relation (labeled by 71) and the 
other is treated as a perturbation (labeled by 72), eventu- 
ally from a different band. It will be considered schemati- 
cally that: 
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Where 71 and 72 keep track of each of the different parts 
of the hopping. The procedure and idea will be clearer and 
useful when working in the momentum space. We will con- 
sider that these terms are characterized by different ranges 
of momenta of holes ^^(ki) and £( 2 )(k2), associated re- 
spectively to the terms E^ 1 ' and Y,( 2 \ such that ki and k2 
can belong to different parts of the band. With this de- 
composition, the following ansatz for the free hole Green's 
functions can be envisaged: 



(G^) y = (d T - n) Sij5(T) + 7 iS (1) y, 
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upper indices ^ stand for the perturbation due to the cor- 
responding part of the hopping term, separated according 
to expression Since this separation is generic and not 
calculated microscopicaly, for the sake of generality we can 
have ^ t^j* depending on the hopping (and band) struc- 
ture. This procedure can be considered such as to provide 
a measure of the (de)localization of the free holes. For in- 
stance, in a normal conducting phase, we should recover 
72 — » 0, that is used in the usual mean field approximation 
|15|lj . On the other hand, at the half filling limit (and very 
close to it) the hopping parameter would be such that 71 = 
0, suitable for the hopping (loop) expansion as discussed in 
details in Ref. pQ. 

3.1. Delocalization of free holes and spin effective action 

We take a continuum limit of the full action, given by ([3]), 
to derive a spin-effective action, S e f / , by integrating out the 
fermion variables with the prescription ([5]) in the same lines 
as it was done in Ref. pQ . Using a finite difference method for 
the term labeled by 71 , we take: £j=j+i — * £(i) + aV£(i) and 
perform a Fourier transformation. Therefore we consider 
free holes are endowed with a dispersion relation e(k( 1 )), 
being that ki (k2)refers to the momenta of modes labeled 
by 71(72). This yields the momentum dependent Green's 
function: Go[/t; e(kW)]. The particular dispersion relation 
e(kW) is completely defined by the lattice (geometry and 
dimensionality). For the sake of the main argument, we con- 
sider a two dimensional square lattice, for which it follows: 
e(k ( i)) ~ 2j(i)tJ2k,a 0k,<7 ( cos ( k x) +cos(k y )) , from what 
the continuum limit is extracted. Prescription ((4| might 
also be associated to a superposition of (nearly) localized 
and (fully) delocalized states. In order simplify the notation 
km momenta will be denoted simply by k from here on. 

The corresponding effective action, with the integration 
of fermions variables, can be written as pQ: 

S eff = TrLogG" 1 = TrLog (G^ 1 - ia + £ (2) ) 

= TYLogGo 1 + TrLog(l - G ia + G S (2) ). (6) 

The different modes of the fermions are decoupled such 
that the corresponding £W are treated (nearly) indepen- 
dently. The free hole Green's function can be calculated, 



for the sake of generality, for a given (sub)lattice A, instead 
of an unique lattice, being written as: (Gq 1 ) a (h a , k) = 
(d T — \i A + eA(k)) _1 (5(r). This case of (at least) two sub- 
lattices will not be worked out here, and this might be con- 
sidered when there are different hoppings in the sublattices 
or between each of them. This can yield the terms labeled 
by 71 and 72 contributing in each of the different sublat- 
tice. The long- wavelength expansion is done by considering 
G 1 >> Z<- 2 \ The reliability of this expansion depends on 
several parameters, seen in expressions (|4|) and ([5]). Basi- 
cally it is required that /i + e(k) >> X^ 2 ) where X^ 2 ) is only 
part of the full hopping term. We remind further that X^ 2 ) 
is basically proportional to t and the long- wavelength limit 
corresponds to a gradient expansion of t < Zi\zj >. There- 
fore we expect to provide a complementary investigation 
to the loop expansion analyzed in Ref . pQ . considering the 
role of the derealization of free holes. For the sake of the 
argument and to show preliminary analytical results, X^ 1 ) 
and X^ 2 ) are considered to somehow decouple from each 
other. In this case the expression for the leading terms of 
the effective action, is obtained in the very same way as 
shown in Ref. [Tj . Keeping track of the time ordering in the 
path integral with a Taylor expansion in fluctuating times 
[181 1] , for D-dimensions S e f / is given by: 
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where ... stands for the higher order terms. In the case the 
spinon dynamics decouples completely from the holes, we 
rewrite S e ff in the momentum space, with zero momen- 
tum transfer between holes and spinons. The effective co- 
efficients J e f f and K e f f can be written as: 
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The K e f f is the coefficient of the time dependent term and 
it will not be analyzed here, since it does not modify the 
magnetic ordering in a first analysis. Kq and K\ are the 
(upper and lower) limiting values of the momenta of holes 
which contribute in that part of the decomposition of the 
variables, i.e. holes from 71 term. For 71 = 0, the expres- 
sions obtained in Ref. [T] are reproduced. In this expression 
the second order (leading) term can be written in the form 
of a Heisenberg coupling (either in k-space or in the lattice) 
as: 



H. 



J, 



eff 



eff 



\(Zi\Zj)? = JeffY,& S i + \) 



where the classical symbols of the spin operators for the 
quantum s = 1/2 Heisenberg model were used |T], and 
the corresponding modes ki and k2 have been separated. 
Below we show some expressions for the coefficient of the 
Heisenberg spin-coupling Jeff- The second order term in 
the effective action (|7|) has different signs and structures 
for J e ff depending on the range of the parameters. How- 
ever it is worth emphasizing that due to this separation, 
the spinon connection (< Zi\zj >) in expression ([7|) corre- 
sponds to only part of the full spinon dynamics (that from 
X^ 2 )), and therefore the resulting phase might not have a 
fully saturated (ferro) magnetic ordering. The separation of 
the full hopping in ranges of momenta, labeled by 71, 72, is 
a relevant assumption for this analysis. A cutoff in the mo- 
mentum integration might correspond to the emergence of 
a kind of Fermi surface for the vacancies in which an inte- 
gration in k is limited by hp, the momentum of holes at the 
(eventual) Fermi surface. A decomposition in low and high 
energy modes suggests a renormalization group analysis, 
which will be presented elsewhere. However separation of 
modes has been also implemented in, for example, Ref. [19]. 

3.2. Some analytical expressions for J e f / 

The particular dispersion relation for the holes was not 
explicitly written so far. For the sake of generality, two cases 
/are considered of the following usual forms: 
efn(k) = 7i5ife = a\k and e(m(k) = 7161k 2 = 61k 2 . 

Changing variables for each of the cases, we write: 
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where ttrj is the integral of the D dimensional solid angle; 
and the cutoffs are: = x[ I] - iii^iM 
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eventual formation of a Fermi surface for the (spinless) va- 
cancies in a normal metallic phase, we could identify the 
chemical potential to: fi = e(kp), where kp is the momen- 
tum at the Fermi surface. The result for the quadratic e(jj) 
in 2-dim is the same as for em in 1-dim, apart from a nor- 
malization. 

The integrations, in 2 dim, yield respectively: 
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A short example can be taken, for T = 0, by choosing 
K\ > Kq ~ kp ~ /i/ai for the first of these expressions, in 
which case iCo might be the momentum at a Fermi surface 
for the holes. The resulting spin-effective coupling at zero 
temperature is given by: J^X ~ 7§i 2 (2ai-RTi — /i) / (Awai) > 
0. This is a ferromagnetic coupling representing a Nagaoka- 
type phase, with some stability since it remains finite for a 
finite corresponding concentration of holes Sh > 0. In par- 
ticular, for very small 5h, we have the compressibility of 
holes, Kh(T — > 0) — > 0. It is interesting to notice that 6h and 
Kh are calculated analytically. As temperature increases the 
decay to a paramagnetic phase should take place, depend- 
ing strongly on the spinon dynamics whose investigation is 
outside the scope of the present work. We emphasize that 
the present work only aim to provide a different starting 
point for investigating the role of delocalization of holes. 
The appearance of the ferromagnetic coupling however was 
related to the range of integration of the momenta of holes. 
Should we consider a different physical picture in which the 
relation among the variables Kq , K\ , \i were related differ- 
ently, it can give rise to an (itinerant) anti-ferromagnetic 
coupling. This analysis remains valid for U < oo. 

In the limit of no limitation on the momenta of holes, i.e. 
Kq = and K\ — > oo, in 2 dim, it yields respectively: 
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We notice that these couplings might depend on fj,. They 
provide anti- ferromagnetic spin Heisenberg couplings. 

4. Final remarks 

We have shown that the delocalization of free holes might 
be a relevant issue for the understanding of the stability 
of the Nagaoka's phase with respect to a finite concen- 
tration of holes. More generally we proposed a framework 
for investigating different magnetic orderings in the limit 
of very large Coulomb repulsion and low concentration of 
holes for the Hubbard model. A spin effective action was 
found to have a form of localized Heisenberg coupling in 
the long wavelength limit along with the work presented 
in Ref. pQ. It can be ferromagnetic (Nagaoka-type phase) 
or anti-ferromagnetic depending on the relation among the 
chemical potential and the eventual values of the limita- 
tion on the summation/integration of momenta carried by 
the holes. For that, the hopping term was separated in two 
parts, corresponding to high and low momentum modes or 
to hopping among different bands 0- A microscopic deriva- 
tion of the prescriptions adopted was not yet presented and 



1 Eventually the parameters 74 introduced to label £00 might be 
expected to depend on the temperature, U and concentration of holes 
&h |20I21I becoming phcnomcnological. For example, one might want 
to account the variation of 7^ due to a finite value for the Coulomb 



for the sake of the main argument it was considered that the 
holes are reasonably decoupled from the spinons. In par- 
ticular, by endowing holes with a dispersion relation such 
that a kind of Fermi surface can be formed, it was found 
that the Nagaoka's phase at finite concentration of holes 
can have some stability in a long- wavelength limit. A more 
complete analysis will be presented elsewhere. 
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repulsion by means of prescriptions. They might be given by: 
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where oti (i = 0, U, 1, 2) depend on the parameters of the model and 
their values are bounded by and 1. The values of such parameters 
must be constrained due to expression l(4). For example for this 
second parameterization, 02 + 03 = 2ajy. This makes possible to 
consider that the hopping term contributes both in S and in Go , 
whereas in the first prescription (i) we recover the development of 
Ref. [I], for which 71 = when U = 00. 
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